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I study the effectiveness of fault-tolerant quantum computation against correlated Hamil- 
tonian noise, and derive a sufficient condition for scalability. Arbitrarily long quantum 
computations can be executed reliably provided that noise terms acting collectively on k 
system qubits are sufficiently weak, and decay sufficiently rapidly with increasing k and 
with increasing spatial separation of the qubits. 
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1 Introduction 

Our planet is in the midst of a digital revolution, validating the scalability of classical infor- 
mation processing. Will quantum computers likewise be scalable, eventually performing tasks 
that surpass what could be done if the world were classical? 

The accuracy threshold theorem for quantum computation establishes that scalability is 
achievable provided that the currently accepted principles of quantum physics hold and that 
the noise afflicting a quantum computer is neither too strong nor too strongly correlated 
[TJ [21 |31 [31 [51 H3 \7\. For scalability to fail as a matter of principle then, either quantum 
mechanics must fail for complex highly entangled systems (as 't Hooft [5] has suggested), or 
else either the Hamiltonian or the quantum state of the world must impose noise correlations 
that overwhelm fault-tolerant quantum protocols (as Alicki et al. [HI QUI E] and Kalai [T21 
US D31HS1 QI] have suggested). 

Because of the profound implications of large-scale quantum computing for computational 
complexity theory and fundamental physics, skepticism is natural and useful. Debate about 
the feasibility of fault-tolerant quantum computation can sharpen our understanding of the 
issues, and raise the stakes as quantum science and technology continue to advance. But 
skeptics should be pressed for a conception of Nature in which classical computing is feasible 
yet quantum computing is forbidden. 

The theory of fault-tolerant quantum computing closely resembles the corresponding clas- 
sical theory, but there are also important differences. Perhaps most fundamentally, a classical 
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computer can perform reliably even if the information being processed leaks to the environ- 
ment, but a quantum computation will fail unless the processed quantum information remains 
almost perfectly concealed. Indeed, in classical systems expelling heat to the environment is 
essential to ensure controllability, while for quantum systems energy dissipation may induce 
decoherence and hence cause trouble. 

A central lesson of fault-tolerant quantum computing is that this is a false dichotomy 
- energy dissipation is just as crucial for reliable quantum computation as for classical 
computation [T7J. The trick is to expel entropy without exposing the protected coherent 
quantum information to the environment. In principle this can be achieved using quantum 
error-correcting codes [HI [19] and carefully designed fault-tolerant protocols [20] [21]. Another 
important lesson is that small rotation errors in imperfect unitary quantum gates can be 
digitized and hence corrected like the bit flip errors in a dissipative classical system. 

The goal of this paper is to exhibit a class of noise models for which quantum computing 
is provably scalable. We will assume that qubits can be refreshed on demand, i.e., that it 
is possible to prepare a standard initial state of n qubits, approximating the product state 
|0)® n , with small, weakly correlated errors. One might imagine that Nature conspires to block 
the creation of a good approximation to a pure product state, but such a limitation would 
threaten the scalability of classical computation as well, and so does not seem like a promising 
way to make a fundamental distinction between classical and quantum computation. 

More plausibly, the distinction might arise because noise correlations unavoidably obstruct 
the creation of profoundly entangled states of many qubits. We will address this issue by 
studying a class of Hamiltonian models of correlated noise and deriving a sufficient condition 
for scalability within the context of this class of models. Skeptics are invited to explain why 
no quantum engineer could ever build a system with noise meeting this criterion. 

In our models, the noise correlations arise from an (in principle infinite) series of terms 
in the Hamiltonian, where terms acting on k system qubits have an operator norm obeying 
an upper bound that drops sufficiently rapidly with increasing k and with growing spatial 
separation among the qubits. It would be desirable to relax our scalability criterion in various 
ways. In particular, scalability can be proven for noise models in which system qubits couple 
to harmonic oscillator bath variables with unbounded norm, assuming the initial state of the 
bath meets certain conditions (for example, if the bath starts out in a low-temperature Gibbs 
state) [22]. But, so far, a proof of scalability for a system in contact with an oscillator bath 
has been worked out only for Gaussian noise, i.e., for the case where the noise is completely 
characterized by the two-point correlations function of the bath variables. Extending that 
argument to a nonlinear oscillator bath with non-Gaussian correlations remains a technically 
challenging open problem, worthy of further attention. 

Gaussian noise models are often regarded as reasonably realistic in physical settings where 
the system is weakly coupled to many environmental degrees of freedom [23) . But even in the 
Gaussian case, arguments for scalability hold only under assumptions about the frequency 
spectrum of bath fluctuations [2U US] [25] [27] , and some skeptics have criticized these 
assumptions 28 , 29 , 30J . However these critics have not clearly identified any class of Gaussian 
noise models which would allow high-fidelity gates in few qubit systems while disallowing 
large-scale fault-tolerant quantum computing. Hence if their objections carry weight, we 
may encounter a barrier blocking further systematic improvements in quantum gate fidelity 
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in the relatively near future. The question we are trying to address here is not whether 
noise will limit the reliability of small-scale quantum computers but rather whether large 
scale quantum computing might eventually fail, even though small scale quantum computers 
continue to improve. 

Anyway, without further apologies, we will use a Hamiltonian model in which the noise 
strength can be characterized using the operator norm. We go beyond previous work [311 
by investigating the effects of not just few-body correlations in the noise but also sufficiently 
weak many-body correlations. The goal is to get a clearer picture how harmful such noise 
correlations could be. The type of model we study has a notable advantage — we do not 
need to make any assumption about the initial state of the bath to derive useful results. We 
formulate the model and state the main result in Sec. 2, then prove it in Sec. 3. 

2 Noise model and scalability criterion 

The noise model we consider is formulated by specifying a time-dependent Hamiltonian H 
that governs the joint evolution of the system and the bath, which can be expressed as 

H = Hs + H b +Hsb; (1) 

here H$ is the time-dependent Hamiltonian of the system that realizes an ideal quantum 
circuit, Hg is the Hamiltonian of the bath, and Hgg, which describes the coupling of the 
system to the bath, is the origin of the noise. We place no restrictions on the bath Hamiltonian 
Hb- Without any loss of generality, we may expand the system-bath Hamiltonian in the form 

h sb = ^2 H ii,i2,...i k = ^2^i Y H iui 2 ,...i k - ( 2 ) 

fe=l (U,l2, • ••«*,) fc=l Uj*2,-"'* 

(k) 

Here, ^ acts on the k system qubits labeled by the indices ii, z 2 , . . . ik, and also acts 

arbitrarily on the bath; for each k we sum over all ways of choosing k system qubits. We 
use (ii, 12, . . . ik) to denote an unordered set of k qubits; by definition, H-*\ ik is invariant 
under permutations of the k qubits and vanishes if two of the indices coincide. Hence the two 
expressions for HgB i n Eq.Q are equivalent. We will not need to assume anything about the 
initial state of the bath, except that the system qubits can be well enough isolated from the 
bath that we can prepare single-qubit states with reasonable fidelity. 

We use the term location to speak of an operation in a quantum circuit performed in a 
single time step; a location may be a single-qubit or multi-qubit gate, a qubit preparation step, 
a qubit measurement, or the identity operation in the case where a qubit is idle during a time 
step. We model a noisy preparation as an ideal preparation followed by evolution governed by 
H, and a noisy measurement as an ideal measurement preceded by evolution governed by H. 
It is convenient to imagine that all system qubits are prepared at the very beginning of the 
computation and measured at the very end; in that case the noisy computation can be fully 
characterized by a unitary evolution operator U acting jointly on the system and the bath, 
obtained by solving the time-dependent Schrodinger equation for the Hamiltonian Eq.Q. 

Let's briefly explain the physical justification for these assumptions. For fault-tolerant 
computing to work, there must be a mechanism for flushing the entropy introduced by noise; 
typically, entropy is removed from the computer by error-correction gadgets which use a 
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supply of fresh ancilla qubits that are discarded after use. For mathematical convenience, 
we suppose that the initial state of the system includes all of the ancilla qubits that will be 
needed during the full course of the computation. To model the actual situation, in which 
ancilla qubits are prepared as needed just before being used, we also suppose that ancilla 
qubits are perfectly isolated from the bath until "opened" at the onset of the gadget in which 
they participate. Similarly, we suppose that the measurements of all ancilla qubits are delayed 
until the very end of the computation, but that these qubits are "closed" (their coupling to 
the bath is turned off) at the conclusion of the gadget in which they participate. Fault- 
tolerant gadgets sometimes also include quantum gates that are conditioned on the classical 
outcomes of earlier measurements. These conditional gates can be included in our framework; 
operations conditioned on measurements may be replaced by coherent gates, conditioned on 
the state of a "closed" control qubit that will be measured later (see Sec. VIC of [H]). With 
these stipulations, our noise model is equivalent to a more realistic one in which ancilla qubits 
are repeatedly measured, reset, and reused. In this model we take for granted that "pretty 
good" fresh ancillas can be prepared at any time, or equivalently that qubits can be effectively 
erased at any time. A similar assumption would be needed to ensure scalability in an analysis 
of fault-tolerant reversible classical computation [T7] . 

Our goal is to derive from Eq.([2| an expression for the effective noise strength e of the noisy 
computation, which is defined as follows [r?l 131). We envision performing a formal expansion 
of U in powers of the perturbation HgBi to all orders. Consider a particular set T r of r circuit 
locations, and let E(I r ) denote the sum of all terms in the expansion such that every location 
in X r is faulty, i.e., such that at least one of the qubits at that location is struck at least once 
by a term in H SB during the execution of the gate. We say that the noise has effective noise 
strength e if 



for any set I r . The accuracy threshold theorem for quantum computing shows that scalable 
quantum computing is possible if e is less than a positive constant Eq ~ 10 -4 [6} I22j. 
Let us define 



where the maximum is over all system qubits and all times, the sum is over all system qubits, 
and £o is the maximal duration of any location. Then our main result can be stated as follows. 

Theorem 1 (Effective noise strength for correlated Hamiltonian noise) If each quantum 
gate acts on at most m qubits and if 



\\E(l r )\\<e r 



(3) 




(4) 



vr < ha k , 



(5) 



for all k, then 




(6) 
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where 

" (k-l)\f k+l (2a) 



OC 



f* = i, + ■ (7) 

1=0 y ' 



It follows that quantum computing is scalable provided the strength of fc-qubit interactions 
decays sufficiently rapidly with fc (so that the sums in Eq.Q and Eq.Q converge), and also 
decays as the spatial separation of the qubits increases (so that the sum defining fj^ in Eq. Q 
converges). 

If, for example, /& = 1, then 

oo 

9k < X>a)' = (I - 2a)- 1 = C(a), (8) 

1=0 

and hence 



e < 2ma 



(e' 6 " 1 '/ 2 ) « 4.72 ma, (9) 



where the last approximation uses C(a) w 1 for a <C 1, as is the case if e is smaller than the 
threshold value £o ~ 10~ 4 . This observation can be restated as the following corollary: 
Corollary 1 If each quantum gate acts on at most m qubits then the effective noise strength 
can be expressed as 

/ , , , \ C(a) 

e < 2ma M 6 " 1 )/ 2 ) , (10) 

where 

l/fc 

a = max (max >J \\H\^ ix Jl *o | , d'H 



xfmax ^ NS 2 , i3 ,..., i J|to] 



t is i/ie maximal duration of any circuit location, and C(a) = (1 — 2a) . 
If instead 

fk < kl/kP (12) 

where p > 1, then 

< fc! /A (fc-l)!fcP/ fc+; (2a)' \ = k[(^ (k-l)\kP(k + l)\(2a) l \ 

9k - kp (k + i-i)\k\ J kp (k + 1 - iy.k\(k + i)p J 



< 



k\ ^kP~ 1 (2a) l \ kl ( 1 \ fc' s 

kP [L jkTiy^ ) - Tp {i^2a ) = kp c{a) - (13) 



J=0 



For p > 1 the sum over fc in Eq.Q converges, and hence we obtain a finite expression for 
e; therefore, scalable fault-tolerant quantum computation is achievable for sufficiently small 
(nonzero) a. We have obtained: 
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Corollary 2 If each quantum gate acts on at most m qubits then for any p > 1 the effective 
noise strength can be expressed as 



e < 2ma p exp 



V fc=l / 



(14) 



where 



a v = max 
F fc>i 



max W H i[ 



1/fe 



(15) 



to is the maximal duration of any circuit location, and C(a) = (1 — 2a) 
In particular, if p = 2 for example, we find 

e < 2ma I e* /12 J w 4.55 ma, (16) 
again using C(a) « 1 to obtain the numerical expression. 
3 Proof of Theorem 1 

In [31] . scalability was proven for the special case in which only the k = 2 term in the 
Hamiltonian is nonzero. To prove Theorem [I] we generalize the ideas used in |31j . We write 
the system-bath Hamiltonian as 



H 



SB 



H SBa 



(17) 



where a is a shorthand for the indices fc, and i\, 12, ■ ■ ■ ik in Eq.([2]). For the sake of conceptual 
clarity we imagine dividing time into infinitesimal intervals, each of width A, and express the 
time evolution operator for the interval (t, t + A) as 



U(t + A, t) ss e 



-iAH 



e -iAHs e -iAH B 



EH 



SB 



iAH; 



SB,a) 



(18) 



(We have omitted terms higher order in ||if||A; strictly speaking, then, to justify Eq.(18) 



we should regulate the bath Hamiltonian Hg by imposing an upper bound on its norm, 
then choose A small enough so these higher order terms can be safely neglected.) We expand 
U(t + A, t) as a sum of monomials, where for each value of a either Igg or —iAHss.a appears; 
then we obtain the perturbation expansion of the full time evolution operator U over time T 
by stitching together T/A such infinitesimal time evolution operators. 

We will refer to the r specified locations in the set I r as the "marked locations" and to the 
remaining locations as the "unmarked locations." For now, suppose for definiteness that all 
of the marked locations are single-qubit gates. For any term in the perturbation expansion 
contributing to E(I r ) there must be an earliest infinitesimal time interval in each of the r 
marked locations where a term H$B,a acts nontrivially on that qubit. Suppose we fix the 
infinitesimal time intervals where these earliest "insertions" of H$b occur, and also fix the 
terms {H$B,a} in the system-bath Hamiltonian that act there, but sum over all the terms 
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in the perturbation expansion acting in other time intervals and on other qubits. Then in 
between the fixed earliest insertions in the marked locations, the joint evolution of the system 
and the bath is governed by a modified Hamiltonian 

(modified) 

^(modified) = Hs + Hb+H {modified) ^ ff (modified) = £ (ig) 

a 

where the modified sum excludes any term Hsb,o, acting nontrivially on any one of the marked 
locations during any time interval prior to the fixed time of the earliest insertion. The im- 
portant point is that the time evolution operator in between successive insertions of the 
perturbation is unitary and hence has unit operator norm. Using the submultiplicative prop- 
erty of the operator norm, then, we conclude that the contribution to E(I r ) with the earliest 
insertions at the marked locations fixed has operator norm bounded above by 

(earliest) 

II (II^b, q ||A), (20) 

a 

where the product is over the terms in the system-bath Hamiltonian that act at the earliest 
insertions. To bound || £7(2^)11, we sum over the to/ A time intervals at each location where 
the earliest insertion may occur, and also sum over all the ways of choosing the term Hsb,u 
that acts at each insertion, obtaining 

(insertions) (earliest) 

\\E(l r \\< ]T J] (W H SB,a\\t ). (21) 

{H S B, a } a 

Summing over the possible intervals for the first insertion turns the factor A into the factor 
t . 

Now we have to figure out how to sum over all ways of choosing the terms {HsB,a} acting 
at the earliest insertions inside the r marked circuit locations. Since Hsb contains multi-qubit 
terms, a single term in H$b can simultaneously produce the first insertion at multiple circuit 
locations occurring in the same time step. Specifically, a single term in might cause 

simultaneous faults in j of the r marked locations for any j < k, if all of these j locations 
occur in the same time step. We use the term "j-contraction" to refer to the case where a 
single term in Eq.Q produces the first insertion in each of j marked locations. See Fig. 1. 

First we find an upper bound on the strength of a one-contraction, the operator norm of 
the sum of all terms that cause one particular circuit location to be faulty. If the qubit at the 
marked location carries the label i\, the term in the Hamiltonian responsible for the earliest 
insertion at this location could be any H^~j^ 3 - 3 Jx for I > 0; here for each m = 1,2,...,/ 
either qubit j m is unmarked or else qubit j m is marked but has already been struck by an 
earlier insertion during the same time step. Hence an upper bound on the strength of the 
one-contraction is 

CO 

m = T,^ 1+l \ (22) 

1=0 
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Fig. 1. Contractions occurring during one computational time step, with time flowing vertically 
downward, where rectangles represent qubits and k dots connected by a horizontal line signify 
an insertion of a term in H^ k \ which acts on k qubits. All six of the qubits shown have faults 
during this time step, arising from a one-contraction, a three-contraction, and a two-contraction. 
The two-contraction is due to a term in that actually afflicts all six qubits, but counts as a 

two-contraction because four of these qubits have already been hit by other contractions earlier in 
the same time step. 



where 

^ i+i H- e ii<a 2 ,.,jMo=4^i e \\H\i2 h ,... t jt, (23) 

Ul,j2,---,jl) jl,j2, — ,jl 

Here, to obtain an upper bound, we sum each index j m over all qubits, whether marked or 
unmarked, and also obtain the factor to by summing over all the infinitesimal time intervals 
during a single time step. We have also maximized this expression over all possible ways to 
choose qubit i\ from among the marked qubits — the superscript "in" in max'" 1 ' indicates 
that we maximize over only the marked qubits, while the superscript "all" in ^( all - ) indicates 
that we sum over all qubits without any restriction. Of course, our upper point would still be 
valid were we to relax the restriction and maximize over all qubits, whether marked or not. 
Similarly, for k > 1, the strength of a ^-contraction can be bounded by 

oo 

m = Y,^ +l \ (24) 

1=0 

where 

^ +0 = E E H«., fcJl ,i 2 ,..,JMo. (25) 

(ii,i 2 ,...,ifc) (ji,j2,---,ji) 

Here, for the "in" sum the qubits are restricted to the marked locations and for the "all" sum 
they may be at either marked or unmarked locations. Note that the upper bound rji involves 
a maximum over marked qubits, while the upper bound rjk for k > 1 involves instead a sum 
over marked qubits; the reason for this distinction is explained in the next paragraph. 

By summing all ways of choosing the first insertion in each of r marked locations, we 
obtain the bound 

(r) oo 

e n^w- ^ 

rx,r-2,T3,... k=l 
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Here is the number of /c-contractions, and the sum J^ 1 *" 1 is subject to the constraint 
Sfc k r k = r - To obtain Eq.(26), we observe that, for k > 1, 

(in) (all) \'' fc 

E E ,J *° ( 27 ) 

K (ii,i 2 ,--,ik) (ji,32,---,jl) / 

contains each way of choosing fc-contractions among the r marked locations r^! times, plus 
additional nonnegative terms; the factor l/r^l in Eq.(26l compensates for this overcounting. 
Furthermore, once all the higher rank contractions have been fixed, the locations where one- 
contractions occur are completely determined. That is why we defined r\\ by maximizing over 
ii, rather than summing i\ over all the marked qubits. 

In Eq.(22) we have derived a bound on the strength of the noise acting at a single circuit 
location. We wish to go further and investigate whether the correlations in noise acting 
collectively on many circuit locations could overcome fault-tolerant protocols, even if the 
individual gates perform very well. For this purpose, we should relate rj^ for k > 1 to 771. 
Note that in for k > 1, we can replace the sum over ways to choose k qubits by a sum 

over all qubits divided by A:!, and similarly we can replace the sum over the ways to choose / 
qubits by a sum over all qubits divided by /!, obtaining 

^ +i) = M. E /j . E \\ H iSZ-i,,h,h,...,J *°- (28) 



By summing i\ over the r marked locations we obtain the bound 

r^fii E \ E NS., fe , 1 , J W ! ll*o; (29) 

»2,— ,ik jl,32,—,jl 

note that we still have a bound if we extend the "in" sum to a sum over all qubits. From 
Eq. ( 23 ) we have 

^Hak^w E E W H iHti h j lda ,...J *o- (30) 
which implies, for k > 1, 

Vt" < - (fc ^f 7 ^r ) = l + \~ l }r i) < r -2^rj{ k+l) - (3D 
Hence we find, for k > 1, 

oo 

^<Y k T. 2k+l ^ +l) - (32) 
1=0 

This is the key inequality that we needed, relating (an upper bound on) the strength of 
collective noise acting on k circuit locations to a sum over (upper bounds on) contributions 
to the noise strength for a single location. 
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Now suppose, as in the hypothesis of Theorem 1, that 

1 *. „k 



(k) 

m = 



(k-iy 



~(k) , fkot 

m < 



(A -1)1 



From Eq. (32 1 we obtain 



where 



5fe 



/* + £ 



(fc-l)!/ fc+ ,(2a)' 
(fc + Z-1)! 



(33) 



(34) 



(35) 



Then the bound Eq.(26) becomes 

1 (rg k (2a) k ^ n 



* r < e n 



2k\ 



-w e ni(i) 



(36) 



ri,f2,)"3,... &=1 v ' r 1 ,r 2 ,r 3 ,... k=l 

recalling the constraint on the sum. If we now relax the constraint on the sum, we have 



oo oo 



<- (M'liE £(!)"-<*■>' IMS) 



k=l r k =0 



k=l 



(2aY exp [J2 



\k=l 



9k 
2fc! 



2a exp ^2 



\k=l 



9k 
2k\ 



(37) 



Up until know, we have considered all circuit locations to be single-qubit locations. In 
the case of an m-qubit gate location, the location is faulty if the system-bath perturbation 
acts nontrivially on any one of m qubits, which enhances each rjk appearing in Eq.(26) by at 
most a factor of m k , and hence increases our upper bound on the effective noise strength by 
at most a factor of m. This completes the proof of Theorem 1. □ 



4 Conclusions and outlook 

Theorem 1, combined with results from [B], shows that fault-tolerant quantum computing is 
scalable in principle for a class of correlated noise models. The key feature of these models is 
that, while there are terms in the noise Hamiltonian acting collectively on k system qubits for 
k S> 1 (and simultaneously on the quantum computer's environment, i.e., the "bath"), the 
operator norm of these terms decays faster than any power of k. The result may apply even 
if, for each fixed fc, the operator norm of the fc-qubit noise term decays algebraically rather 
than exponentially as the qubits are spatially separated, provided the decay is sufficiently 
rapid as a function of distance for the sum over system qubits in Eq.Q to converge. 

Theorem 1 generalizes results found in [31] [5S] for the case k = 2. Though the analysis 
is formulated in terms of interaction-picture perturbation theory, it is rigorous because our 
estimate of the effective noise strength is derived by bounding perturbation theory summed 
to all orders. 
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We also assume that at all times during the computation it is possible to prepare qubits 
in a standard initial state and to measure qubits in a standard basis with reasonable fidelity, 
but we make no other assumptions about the state of the bath. Although the bath might 
be quite "hot," the flow of entropy from the bath to the computer is impeded by the weak 
coupling between the system and the bath, allowing a fault-tolerant protocol to maintain a 
steady state where entropy is removed fast enough to keep the computer "cool." Entropy is 
carried away by the ancilla qubits used in error correction gadgets, which are subsequently 
erased and reused. A mechanism for flushing entropy is required in any scheme for stabilizing 
a noisy computer, whether quantum or classical 17 \. 

Of course, the condition for scalability derived in Theorem 1 is merely sufficient and not 
necessary. In particular, scalability may be provable even if the system-bath Hamiltonian 
has unbounded norm, but in that case further assumptions are needed about the state of the 
bath at the beginning of the computation. For example, threshold theorems for Gaussian 
noise were proven in [22] , which apply if the bath is a linear system of harmonic oscillators 
and the initial state of the bath is Gaussian (a thermal state for example). In that case, 
the criterion for scalability can be stated as a property of the bath's two-point correlation 
function, ensuring that spatial and temporal noise correlations decay sufficiently rapidly. In 
this particular setting, at least, we can address the concern expressed in [9] that over the 
course of a long computation an oscillator bath could be driven to a highly adversarial state 
that overwhelms fault-tolerant protocols. 

The purpose of this work is to address an issue of principle concerning the scalability of 
quantum computers subject to correlated noise, not to obtain optimized realistic estimates of 
the accuracy threshold for quantum computing. Indeed the noise strength e appearing in our 
criterion for scalability is in effect an error amplitude per gate rather than an error probability 
per gate, and must be below about 10 -4 for known threshold theorems to apply [51 [55]. This 
criterion is probably much too pessimistic — for the general class of Hamiltonian noise models 
considered here, we cannot easily rule out substantial enhancement of the logical failure 
probability due to constructive interference of many coherently combined fault histories, even 
though this seems quite unlikely in practice. Furthermore, we have not attempted here to 
assess the effectiveness of methods such as noiseless subsystems [32] or dynamical decoupling 
[33] which could suppress the noise correlations. Characterizing the residual noise correlations 
when dynamical decoupling is employed seems difficult for general Hamiltonian noise models, 
though some preliminary steps were reported in |34j . Also, to derive Theorem 1 we made 
no assumptions about the bath Hamiltonian Hb, and it may be possible to derive stronger 
results contingent on physically motivated limitations on the bath dynamics such as locality 
constraints. 

The modest results derived here can hardly be expected to assuage the quantum comput- 
ing skeptics, but may nevertheless help to clarify the debate. Can we identify fundamental 
principles of physics that are compatible with large-scale classical computing yet incompatible 
with large-scale quantum computing? Enlarging the class of noise models for which quantum 
computing is provably scalable should shed light on this fascinating and important question. 
If quantum mechanics is valid and if future quantum engineers can devise a controllable many- 
qubit system with noise meeting the criterion derived in this paper, then reliable large-scale 
quantum computing will be achievable. 
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